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Abstract 

Equation of State (EoS) for hot and dense nuclear matter with a quark-hadron 



O ■ phase transition is constructed within a statistical mixed-phase model assuming 

coexistence of unbound quarks in nuclear surrounding. This model predicts the 
Zj ', deconfinement phase transition of the crossover type. The so-called "softest point" 

O j effect of EoS is analyzed and confronted to that for other equations of state which 

CT) • exhibit the first order phase transition (two-phase bag model) or no transition at 

all (hadron resonance gas). The collective motion of nucleons from high-energy 
heavy-ion collisions is considered within a relativistic two-fluid hydrodynamics for 
k>" \ different EoS. It is demonstrated that the beam energy dependence of the directed 

^ | flow is a smooth function in the whole range from SIS till SPS energies and allows to 

disentangle different EoS, being in good agreement with experimental data for the 
statistical mixed-phase model. In contrast, excitation functions for relative strange- 
rs ■ ness abundance turn out to be insensitive to the order of phase transition. 



Key-words: QCD phase transition, heavy-ion collisions, hydrodynamics, directed 
flow, strange particle production. 



1 Introduction 

The predicted phase transition from confined hadrons to a deconfined phase of their con- 
stituents (i.e. the asymptotically-free quarks and gluons or the so-called Quark-Gluon 
Plasma, QGP) is a challenge to the theory of strong interaction. Over the past two decades 
a lot of efforts has been spent to both the theoretical study of deconfinement phase transition 
and the search for its possible manifestation in relativistic heavy ion collisions, properties 
of neutron stars and Universe evolution. A unique opportunity provided by relativistic 
heavy-ion collisions allowed to reach a state with temperature and energy density exceed- 
ing the critical values, T c ~ 170 MeV and e c ~ 1 GeV/ fm 3 , specific for the deconfinement 
phase transition. A rather long list of various signals for the QGP formation in hot and 
dense nuclear matter is available now and it has been probed in experiments with heavy 
ions. Unfortunately, there is no crucial signal for unambiguous identification of the decon- 
finement phase and, for a particular reaction at the given bombarding energy, practically 
every proposed signal can be simulated to some extent by hadronic interactions. 

In this paper we turn to the study of excitation functions for observables to be sensitive 
to the expected QCD deconfinement phase transition. Its manifestation has been considered 
already some time ago by fll], [| . Since a phase transition slows down the time evolution of 
the system due to softening of the EoS, and one can expect a remarkable loss of correlations 
around some critical incident energy resulting in definite observable effects. 
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2 Equation of state in mixed phase model 



Following the common strategy of the two-phase (2P) bag model ||, one can determine 
the deconfinement phase transition by means of the Gibbs conditions matching the EoS 
of a relativistic gas of hadrons and resonances, whose interactions are simulated by the 
Van der Waals excluded volume correction, to that of an ideal gas of quarks and gluons, 
where the change in vacuum energy in a QGP state is parameterized by the bag constant 
B. Thermodynamics of the hadron gas is described in the grand canonical ensemble. All 
hadrons with the mass rrij < 1.6 GeV i are taken into consideration. One should emphasize 
that the phase transition in the 2P model is right along of the first order by constructing. 

To reproduce the variety of phase transitions predicted by QCD lattice calculations we 
represent a phenomenological Mixed Phase (MP) model [^j, The underlying assumption 
of the MP model is that unbound quarks and gluons may coexist with hadrons forming a 
homogeneous quark/gluon-hadron phase. Since the mean distance between hadrons and 
quarks/gluons in this mixed phase may be of the same order as that between hadrons, the 
interaction between all these constituents (unbound quarks/gluons and hadrons) plays an 
important role and defines the order of the phase transition. 

Within the MP model (§ g the effect ive Hamiltonian is expressed in the quasipar- 
ticle approximation with density-dependent mean-field interactions. Under quite general 
requirements of confinement for color charges, the mean-field potential of quarks and gluons 
is approximated by 

U q (p) = U g (p) = ^ ; 7 >0 (1) 
with the total density of quarks and gluons 

P = Pi + Pg + Y, v 3Pj > 

3 

where p q and p g are the densities of unbound quarks and gluons outside of hadrons, while pj 
is the density of hadron type j and z/j is the number of valence quarks inside. The presence 
of the total density p in (|l|) implies interactions between all components of the mixed phase. 
The approximation ([]]) mirrors two important limits of the QCD interaction. For p — > 0, 
the interaction potential approaches infinity, i.e. an infinite energy is necessary to create 
an isolated quark or gluon, which simulates the confinement of color objects. In the other 
extreme case of large energy density corresponding to p — > oo, we have U q = U g — which 
is consistent with asymptotic freedom. 

The use of the density-dependent potential ([!]) for quarks and the hadronic potential, 
described by a modified non-linear mean- field model ||, requires certain constraints to be 
fulfilled, which are related to thermodynamic consistency [|], |^]. For the chosen form of the 
Hamiltonian these conditions require that U g (p) and U q (p) do not depend on temperature. 
From these conditions one also obtains a form for the quark-hadron potential M. 

A detailed study of the pure gluonic £77(3) case with a first-order phase transition allows 
one to fix the values of the parameters as 7 = 0.92 and A 1/{37+1) = 250 MeV. These values 
are then used for the £77(3) system including quarks. As is shown in Fig.l for the case of 
quarks of two light flavors at zero baryon density (% = 0), the MP model is consistent 
with lattice QCD data providing a continuous phase transition if the crossover type with 
a deconfinement temperature Td ec — 153 MeV. For a two-phase approach based on the 
bag model a first-order deconfinement phase transition occurs with a sharp jump in energy 
density e at Tdec close to the value obtained from lattice QCD. 
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Figure 1: The reduced energy density and pressure (the Esb an d Psb are corresponding 
Stephan-Boltzmann quantities) of the 577(3) system with two light flavors for ub = 
calculated within the MP (solid lines) and bag (dashed lines) models. Circles and squares 
are lattice QCD data obtained within the Wilson and Kogut-Susskind [Q schemes. 
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Figure 2: The (e,p/e)-representation of the EoS for the two-flavor SU(3) system at various 
baryon densities %. Notation of data points and lines is the same as in Fig.l. 



Though at a glimpse the temperature dependencies of the energy density e and pressure 
p for the different approaches presented in Fig.l look quite similar, there is a large difference 
revealed when p/e is plotted versus e (cf. Fig. 2, left panel). The lattice QCD data differ 
at low e, which is due to difficulties within the Kogut-Susskind scheme || in treating the 
hadronic sector. A particular feature in the MP model is that, for % = 0, the softest 
point of the EoS, defined as a minimum of the function p(e)/e U, is not very pronounced 
and located at comparatively low values of the energy density: Esp ~ 0.45 GeV/fm 3 , 
which roughly agrees with the lattice QCD value 0. This value of e is close to the energy 
density inside a nucleon, and hence, reaching this value indicates that we are dealing with a 
single big hadron consisting of deconfined matter. In contradistinction, the bag-model EoS 
exhibits a very pronounced softest point at large energy density Esp ~ 1.5 GeV/fm 3 |9|, p~0| ] . 

The MP model can be extended to baryon-rich systems in a parameter-free way 0, [J. 
As demonstrated in Fig. 2 (right panel), the softest point for baryonic matter is gradually 
washed out with increasing baryon density and vanishes for ns > 0.3 uq (uq is normal 
nuclear matter density). This behavior differs drastically from that of the two-phase bag- 
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model EoS, where esp is only weakly dependent on % |], llOfl . It is of interest to note 
that the interacting hadron gas model has no softest point at all and, in this respect, its 
thermodynamic behavior is close to that of the MP model at high energy densities 0] . 

These differences between the various models of EoS should manifest themselves in 
dynamics discussed below. 



3 Directed flow of baryons 

The EoS described above is applied to a two- fluid (2F) hydrodynamic model |TI| , which 



takes into account finite stopping power of colliding heavy ions. In this dynamical model, 
the total baryonic current and energy-momentum tensor are written as 

r = Jg + J? , (2) 

= r^ + rr , (3) 

where the baryonic current J£ = n a u^ and energy-momentum tensor T£ v of the fluid a 
are initially associated with either target (a = t) or projectile (a = p) nucleons. Later on 
these fluids contain all hadronic and quark-gluon species, depending on the model used 
for describing the fluids. The twelve independent quantities (the baryon densities n a , 4- 
velocities normalized as u ail u^ = 1, as well as temperatures T and pressures p of the 
fluids) are obtained by solving the following set of equations of two-fluid hydrodynamics [JITl 



where the coupling term 



n p n t 



« = , (4) 

W = K , (5) 



(v re i J da NN ^ NX (s) (p-p a y\ (6) 



characterizes friction between the counter-streaming fluids. The cross sections da^N-^NX 
take into account all elastic and inelastic interactions between the constituents of different 
fluids at the invariant collision energy s 1 / 2 with the local relative velocity V re i = [s(s — 
Am 7 N )] 1 / 2 /2m 2 N . The average in (|6|) is taken over all particles in the two fluids which are 
assumed to be in local equilibrium intrinsically |IT| . The set of Eqs. (§) and (|) is closed 
by EoS, which is naturally the same for both colliding fluids. 

Following the original paper PH , it is assumed that a fluid element decouples from the 
hydrodynamic regime, when its baryon density Ub and densities in the eight surrounding 
cells become smaller than a fixed value rif. A value rif = 0.8n is used for this local freeze- 
out density which corresponds to the actual density of the freeze-out fluid element of about 
0.6 - 0.7 n . 

The directed flow characterizes the deflection of emitted hadrons away from the beam 
axis within the reaction x — z plane. In particular, one defines the differential directed flow 
by the mean in-plane component (p x {y)) of the transverse momentum at a given rapidity y. 
This deflection is believed to be quite sensitive to the elasticity or softness of the EoS and 
can be quantified in two ways: In terms of the derivative (a slope parameter) at mid-rapidity 



F _ d (p x (y)) 
dy 



(7) 

y=ycm 
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which is quite suitable for analyzing the flow excitation function, and by another integral 
quantity to be less sensitive to possible rapidity fluctuations of the in-plane momentum: 

d 3 N\ 



J dp x dp y dy p x \E ^ 
J dp x dp y dy \ E ~ dp 



where the integration in the cm. system runs over the rapidity region [0, y cm \- Excitation 



functions in the SIS-AGS-SPS energy range are plotted in Fig.3 for both characteristics [12 
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Figure 3: Excitation functions of the slope parameter F y (left panel) and the average 
directed flow (right panel) for baryons from Au + Au collisions within hydrodynamics 
and different transport simulations. Collected experimental points are taken from [|12]1 . 
The results of transport calculations for three different codes (left lower panel) are given 
by the thin solid (RQMD), dashed (ARC) and dot-dashed (ART) lines (cited according 
to ]13[). The solid line (RBUU) is taken from ]TJ]]. 2F hydrodynamics with the MP EoS 
at the impact parameter 3 fm is compared with the corresponding results of IF- |Kj (right 
upper panel) and 3F- (right lower panel) |Tj| hydrodynamics with the bag-model EoS. IF 
calculations both with and without the phase transition (PT) are displayed. 



Our first 2F hydrodynamic calculations of F y (Ei a b) are in a good agreement with ex- 
periment in the whole energy range considered. In the left lower panel of Fig.3 our results 
are compared with transport calculations. The ARC and ART are cascade models, while 
the RQMD takes also into account mean-field effects. Though all these models agree with 
experimental data at E[ a ^ « 10 A-GeV (considered as a reference point), values of F y at 
lower energies are clearly underestimated, as is evident from comparison with results of 
the E895 Collaboration [l(| (see empty squares in Fig.3). Recently, a good description 
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of experimental points (including the E895 data) was reported within a relativistic BUU 
(RBUU) model [|14]]. The good agreement with experiment was achieved here by a special 
fine tuning of the mean fields involved in the particle propagation. 

The calculated excitation functions of (P x ) for baryons within different hydrodynamic 
models are shown in the right panel of Fig. 3. Conventional IF hydrodynamics for pure 



hadronic matter |10j results in a very large directed flow due to the inherent instantaneous 



stopping of the colliding matter. This instantaneous stopping is unrealistic at high beam 
energies. If the deconfinement phase transition, based on the bag- model EoS fl(| , is included 
into this model, the excitation function of (P x ) exhibits a deep minimum near Ei a t> « 6 
A-GeV, which manifests the softest-point effect of the bag-model EoS as shown in the right 
panel of Fig.2. 

The result of 2F hydrodynamics with the MP EoS noticeably differs from the IF cal- 
culations. After a maximum around 1 A-GeV, the average directed flow decreased slowly 
and smoothly. This difference is caused by two reasons. First, as follows from Fig.2, the 
softest point of the MP EoS is washed out for > 0.4. The second reason is dynamical: 
the finite stopping power and direct pion emission change the evolution pattern. The latter 
point is confirmed by comparison to three-fluid calculations with the bag EoS plotted 
in the right lower panel of Fig. 3. The third pionic fluid in this model is assumed to interact 
only with itself neglecting the interaction with baryonic fluids. Therefore, with regard to 
the baryonic component, this three-fluid hydrodynamics |15, 17] is completely equivalent 
to our two-fluid model and the main difference is due to the different EoS. As seen in Fig. 3, 
the minimum of the directed flow excitation function, predicted by the one-fluid hydro- 
dynamics with the bag- model EoS, survives in the three-fluid (nonunified) regime but its 
value decreases and its position shifts to higher energies. If one applies the unification pro- 
cedure of [|15|], which favors fusion of two fluids into a single one, and thus making stopping 
larger, three-fluid hydrodynamics practically reproduces the one-fluid result and predicts 
in addition a bump at Ei a b « 40 A-GeV. 



4 Strangeness production 

Enhanced strangeness production as compared to proton-proton or proton-nucleus colli- 
sions is one of the QGP signals proposed a long time ago. In the hydrodynamic model 
described above only baryon charge rather than strangeness exchange is included. So, to 
see experimental consequences of EoS with different phase transition order, we consider 
an expanding homogeneous blob of the compressed and heated QCD matter (a fireball) 
formed in heavy-ion collisions. The initial state (eq and ns) for this fireball is estimated 
from results of the QGSM transport calculations in the center-of-mass frame inside a cylin- 
der of the volume Vo with radius R = 5 fm and Lorentz-contracted length L = 2R/j c _ m _. 
Isoentropic expansion is treated in an approximate manner assuming V ~ Vot until the 
freeze-out point defined by Sj = 0.15 GeV/ fm 3 ~ m N n (for more detail see 



One should note that till this point the grand canonical ensemble was used where com- 
plete chemical equilibrium is assumed and the strangeness conservation is controlled on 
average by the strange chemical potential us- In the thermodynamical limit, fluctuations 
in a number of strange particles are small and coincide with those for the canonical ensem- 
ble. However, it is not the case for finite systems at relatively small T where the strangeness 
canonical ensemble should be applied, taking into account the associative nature of strange 
particle creation by exact and local conservation of strangeness. Using the general formal- 
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ism for the canonical strangeness conservation proposed in [19, 20], the partition function 
of a gas of hadrons with strangeness s» = 0, ±1, ±2, ±3 and total strangeness 5 = can be 
written as follows 

1 



— / _d<j> exp(£ S a e»*) 



(9) 



where S s = VJ2i%i • Here is the one-particle partition function for species i and the 
sum is taken over all particles and resonances carrying strangeness Sj. A number of strange 
particles can be found by the appropriate differentiation of the partition function Z$, Eq.(|5]). 
It is easy to see that canonical result can be obtain in the Boltzmann approximation from 
the grand canonical one by replacing the strange fugacity in the following way : 



exp(/i s /T) 



Si 



x 



I (x) 



(10) 



where the argument of the Bessel functions x = 2y / <Si«S_i ~ V. This receipt was applied to 
our treatment of particle abundance at the freeze-out point. Generally, the correlation vol- 
ume in the suppression factor I s (x)/Iq(x) of ( fLOf) does not coincide with the system volume 
V. In our model, the initial Lorentz-contracted volume Vq is considered as a strangeness 
correlation volume. 
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Figure 4: Particle ratios for strange hadrons in full An angle interval for central Au + Au 
collision as a function of bombarding energy. The compilation of available experimental 
points is taken from [f2T], 22|. The calculated excitation functions represent four modeling 
EoS with the canonical suppression factor I s (x)/Iq(x). For the case of the MP model, the 
grand canonical results (dashed lines) are given, as well. 
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Inspection of Fig.|] shows that the inclusion of the canonical strangeness suppression 
factor allows one to decrease noticeably strange particle abundance. However, comparing 
canonical and grand canonical results for the MP model, one can see that they do not 
coincide at high energies as it would be expected. This is explained by the beam energy 
dependence of the strangeness correlation volume in contrast with usual canonical descrip- 
tion [19, 20]. The most striking result followed from Fig.[|is that all the models considered 



in the strangeness canonical ensemble predict practically the same relative abundance of 
strange hadrons in the whole energy range studied. The measured excitation function for 
K + /ir + are reproduced reasonably well excluding maybe the SIS energy. In the case of 
K~ /ir~ the general form of excitation functions also agrees with experimental one but the 
relative abundance is overestimated what mainly originates from neglecting the electric 
charge (isospin) conservation. A simple estimate shows that taking into account the isospin 
conservation the K~ /ir~ ratio decreases by about 18% and 12% at Ei a b = 10 and 150 AGeV, 
respectively, without any essential influence on the K + /tt + ratio. The relative yield of hy- 
perons, in particular A/7r + 's, seems to be overshot in the energy range E\ a \, < 10 AGeV 
what can result from the simplified dynamical treatment: The Bjorken-like longitudinal 
expansion can be applied at the SPS energies, but at the SIS energies the transverse ex- 
pansion is not negligible. It is worthy to note that all calculations have been done with the 
same shock-like freeze-out condition for every EoS without any special tuning. 



5 Conclusions 

It has been shown that the directed flow excitation functions F y and (P x ) for baryons are 
sensitive to the EoS, but this sensitivity is significantly masked by nonequilibrium dynamics 
of nuclear collisions. Nevertheless, the results indicate that the widely used two-phase EoS, 



based on the bag model || |10| and giving rise to a first-order phase transition, seems to be 
inappropriate. The neglect of interactions near the deconfinement temperature results in an 
unrealistically strong softest-point effect within this two-phase EoS. Smooth experimental 
excitation function of the directed flow is reasonably reproduced within the MP model. 

The dynamical trajectories for a fireball state in the T — [is plane are quite different for 
different EoS |18|, However, after using the shock-like freeze-out, the global strangeness 
production is turned out to be completely insensitive to the particular EoS as illustrated 
by the calculation of excitation functions for K + ,K~ and hyperons. To get agreement 
with experiment the canonical suppression factor should be taken into account. The only 
trace of dynamics is the beam-energy dependence of the strangeness correlation volume. 
This effect results in a negative slope of the K + /tt + excitation function at high energies 
which is not reproduced by the equilibrium statistical model with canonical account of 
strangeness |20|, pl[ . 

Among other signals of the QCD deconfinement phase transition, the dilepton and hard 
photon production is the most promising. Being sensitive to the whole evolution time of 
a system, these signals can become experimentally observable to disentangle the crossover 
phase transition, predicted by the MP model, from the first order one which is peculiar for 
the bag-model EoS. 

Useful and numerous discussions with B. Friman, Yu.B. Ivanov, E.G. Nikonov, W. Noren- 
berg and K. Redlich are acknowledged. This work was supported in part by DFG (project 
436 RUS 113/558/0) and RFBR (grant 00-02-04012). 
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